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CYCLIC  DEBONDING  OF  UNIDIRECTIONAL  COMPOSITE  BONDED  TO 
ALUMINUM  SHEET  FOR  CONSTANT -AMPLITUDE  LOADING 


George  L.  Roderick,*  Richard  A.  Everett,  Jr.,* 
and  John  H.  Crews,  Jr. 

Langley  Research  Center 

SUMMARY 

Cyclic  debonding  rates  were  measured  during  constant -amplitude  loading  of  speci¬ 
mens  made  of  graphite/epoxy  bonded  to  aluminum  and  S-glass/epoxy  bonded  to  aluminum. 
Both  room -temperature  and  elevated-temperature  curing  adhesives  were  used.  Debond¬ 
ing  was  monitored  with  a  photoelastic  coating  technique.  The  debonding  rates  were  com¬ 
pared  with  three  expressions  for  strain-energy  release  rate  calculated  in  terms  of  the 
maximum  stress,  stress  range,  or  a  combination  of  the  two. 

The  debonding  rates  were  influenced  by  both  adherend  thickness  and  the  cyclic 
stress  ratio.  For  a  given  value  of  maximum  stress,  lower  stress  ratios  and  thicker 
specimens  produced  faster  debonding.  Microscopic  examination  of  the  debonded  sur¬ 
faces  showed  different  failure  mechanisms  both  for  identical  adherends  bonded  with  dif¬ 
ferent  adhesive  and,  indeed,  even  for  different  adherends  bonded  with  identical  adhesives. 

The  expressions  for  strain -energy  release  rate  correlated  the  data  for  different 
specimen  thicknesses  and  stress  ratios  quite  well  for  each  material  system,  but  the 
form  of  the  best  correlating  expression  varied  among  material  systems.  Consequently, 
empirical  correlating  expressions  applicable  to  one  material  system  may  not  be  appro¬ 
priate  for  another  system,  and  caution  should  be  exercised  in  drawing  parallels  between 
different  material  systems  without  supporting  test  data. 

INTRODUCTION 

Adhesive  bonds  are  widely  used  to  join  structural  components.  Such  joints  produce 
milder  stress  concentrations  than  do  mechanically  fastened  joints.  Adhesively  bonded 
joints  are  particularly  advantageous  for  composite  materials  because  stress  concentra¬ 
tions,  such  as  bolt  holes,  can  significantly  reduce  the  static  strength  of  such  materials. 

Bonds  can  also  be  used  to  join  layers  of  different  materials  to  make  one  new, 
improved  material.  For  example,  hybrid  bonded  systems  of  metal  and  composite 
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layers  are  stronger  for  equal  weight  and  stiffness  than  metals  alone  (ref.  1).  Indeed, 
composites  themselves  derive  their  high  efficiency  from  a  collection  of  bonded  constit¬ 
uent  materials.  However,  under  cyclic  loading  the  bond  is  often  the  weakest  point  in  a 
bonded  system.  Once  started,  cyclic  debonding  usually  continues  progressively  (ref.  2) 
and  compromises  the  integrity  of  the  bonded  system.  Consequently,  when  debonding  is 
detected  in  a  structure,  its  rate  of  propagation  is  crucial  in  determining  the  component 
reliability  and  in  establishing  inspection  intervals.  Unfortunately,  very  little  is  known 
about  cyclic  debonding,  and  a  designer  has  no  rationale  to  predict  debonding  rates.  A 
concept  which  can  be  used  as  the  basis  of  such  a  rationale  is  presented  in  this  paper  and 
is  used  to  analyze  cyclic  debonding  data. 

This  paper  presents  an  analysis  for  constant -amplitude  cyclic  debonding  of  simple 
laminated  specimens  made  of  aluminum  alloy  sheets  bonded  to  a  unidirectional  graphite 
or  fiberglass  composite.  The  analysis  is  based  on  correlations  of  the  observed  debond¬ 
ing  rates  with  calculated  rates  of  strain-energy  release  that  accompanies  cyclic  debond¬ 
ing.  Because  the  debonding  rates  are  intimately  related  to  the  debonding  failure  modes, 
typical  failure  surfaces  are  examined  for  each  of  the  three  material  systems  studied. 

SYMBOLS 

The  units  for  the  physical  quantities  defined  in  this  paper  are  given  in  the  Inter¬ 
national  System  of  Units  (SI)  (ref.  3).  The  measurements  and  calculations  were  made 
in  the  U.S.  Customary  Units. 

a  debond  length,  m 

da/dN  debond-propagation  rate,  m/cycle 

c,n  curve -fit  parameters 

E  YoungTs  modulus,  Pa 

G  strain -energy  release  rate,  j/m 

L  length,  m 

N  number  of  cycles 

P  applied  load,  N 

R  ratio  of  minimum -to -maximum  applied  stress 
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RT 


room  temperature 


r  residual  from  least-squares  curve  fit,  m/cycle 

S  stress  in  region  A  of  composite  sheet,  Pa 

t  thickness,  m 

AT  change  in  temperature  between  cure  and  test  temperatures,  K 

U  strain  energy,  J 

V  volume,  m^ 

w  specimen  width,  m 

x,y  Cartesian  coordinates,  m 

a  thermal  expansion  coefficient,  K~l 

5  deflection,  m 

e  strain 

o*  stress,  Pa 

<j)  strain-energy  density,  j/m^ 

Subscripts: 

A  region  A 

B  region  B 

C  region  C 

max  maximum  value  of  stress 

min  minimum  value  of  stress 
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range  range  of  stress 

max-range  product  of  maximum  stress  and  range  of  stress 

1  aluminum 

2  composite 


EXPERIMENTAL  PROCEDURE 
Specimens  and  Loading 

Figure  1  shows  the  specimen  configuration  used  in  the  present  study.  The  speci¬ 
men  consisted  of  two  aluminum  alloy  sheets  adhesively  bonded  to  a  unidirectional  com¬ 
posite  core.  (See  table  I  for  the  material  properties.)  Three  combinations  of  materials 
were  tested  and  will  be  referred  to  as  material  systems  I,  II,  and  III,  as  shown  in  the 
following  table: 


Material 

system 

Composite 

Adhesive 

Type 

Thickness, 

mm 

Cure 

temperature, 

K 

i 

Graphite/epoxy 

EA-927Ra 

0.13 

294  (RT) 

H 

Graphite/epoxy 

AF-126b 

.13 

394 

III 

S -glass/epoxy 

AF-126b 

.13 

394 

a In  EA-927R,  the  scrim  cloth  is  a  woven  cloth  with  fibers  parallel  and  per¬ 
pendicular  to  the  longitudinal  axis  of  the  specimen. 

h  in  AF-126  the  scrim  cloth  is  a  mat  cloth  with  randomly  oriented  fibers. 


The  specimen  was  designed  for  rapid  debond  initiation.  A  severe  stress  concentration 
was  introduced  in  the  form  of  the  abrupt  change  of  section.  Under  cyclic  loading,  debond¬ 
ing  started  readily  at  this  stress  concentration. 

The  experimental  program  included  tests  for  ranges  of  material  thicknesses,  maxi¬ 
mum  applied  stresses,  and  stress  ratios.  (See  tables  II,  III,  and  IV.)  For  all  the  tests, 
the  cyclic  frequency  was  10  Hz. 

Measurement  of  Debonding  Rates 

A  photoelastic  technique  was  used  to  monitor  the  debond  front.  Figure  2  shows  the 
location  of  the  photoelastic  coatings  on  the  specimen  and  typical  isochromatic  fringes. 
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Photoelastic  coatings  were  bonded  to  the  aluminum  surfaces  of  the  specimen;  the  sur¬ 
faces  were  viewed  through  a  polarizer  and  quarter -wave  plate.  Under  load,  isochromatic 
fringes  developed  at  a  debond  front  as  a  result  of  the  high  strain  gradient  in  that  vicinity. 
The  test  specimens  were  photographed  at  specified  intervals  throughout  each  test  to 
record  the  position  of  the  debond  front. 

STRAIN-ENERGY  RELEASE  RATE  EQUATIONS 


Basic  Equation 

Fatigue -crack-propagation  rates  in  metals  have  been  correlated  with  the  strain- 
energy  release  rate  G  (ref.  4),  where 


G  =  P  —  -  — 
da  da 


(i) 


The  term  P(d5/da)  is  the  work  done  by  the  applied  load  as  the  crack  extends,  and 
dU/da  is  the  change  in  the  stored  strain  energy  as  the  crack  extends.  Conceptually, 
because  debond  propagation  in  composites  is  analogous  to  fatigue -crack  propagation  in 
metals,  such  a  correlation  may  also  be  valid  for  debond  propagation.  Accordingly,  an 
expression  for  G  was  developed  for  the  specimen  configuration  tested.  This  expres¬ 
sion  was  determined  using  a  one -dimensional  elasticity  analysis  described  in  the  follow¬ 
ing  discussion. 

To  simplify  the  analysis,  the  specimen  was  divided  into  three  regions  of  lengths  L^, 
Lg,  and  (fig.  3).  In  regions  A  and  C,  only  uniform  strains  in  the  x-direction  were 

significant.  Consequently,  these  regions  could  be  analyzed  by  an  elementary  elasticity 
method.  In  contrast,  region  B  had  a  complex  strain  distribution  and  could  not  be  ana¬ 
lyzed  using  elementary  methods.  But  the  strain  distribution  in  region  B,  and  conse¬ 
quently  the  length  of  region  B,  can  be  expected  to  remain  constant  as  the  debond  extends. 

As  will  be  shown,  this  constancy  eliminated  the  need  to  calculate  the  strain  distribution 
in  region  B. 

An  expression  for  d<5/da  in  equation  (1)  was  derived  from  the  increase  in  the  end 
deflection  of  the  specimen  dS  caused  by  an  increment  of  debonding  da.  Before  an 
increment  of  debonding  occurred,  the  end  deflection  was  given  by 


6  -  5A  +  5B  +  6C  (2) 

and  after  debonding,  by 

5'  =  8'a  +  5b  +  8'c  (3) 
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Because  the  strain  distribution  in  region  B  was  assumed  to  be  the  same  before  and  after 
the  region  translated,  the  term  did  not  change.  To  find  A5,  equation  (2)  was  sub¬ 

tracted  from  equation  (3);  thus, 

6’  -  6  =  A5  =  6^  -  5A  +  6^  -  6C  (4] 

The  deflections  on  the  right  side  of  equation  (4)  can  be  expressed  in  the  general  form 

6=  eh  (5; 

Substitution  of  equation  (5)  into  equation  (4)  yields: 


A5  —  €a(La  +  Aa)  -  €aLa  +  €/~i (Lo  -  Aa)  -  60 L, 


A\  A  ^  -  A  A  TCCrC 


A6=  eA  -  ec  Aa 


lim  _e  (8 

Aa— 0  Aa  da  A  C 

The  expression  for  dU/da  in  equation  (1)  was  derived  by  calculating  the  change 
of  strain  energy  dU  in  the  specimen  resulting  from  an  increment  of  debonding  da. 
Employing  the  same  reasoning  used  in  the  development  of  equation  (4)  allows  the  change 
in  strain  energy  to  be  expressed  as 

UT  -  U  =  AU  =  UA  -  UA  +  -  Uc  (9 

The  strain  energies  on  the  right  side  of  equation  (9)  can  be  expressed  in  the  general 
form 

U  =  </>V  (10 


where  <p  is  the  strain-energy  density  and  V  is  the  volume  of  strained  material.  Sub 
stitution  of  equation  (10)  into  equation  (9)  yields 

AU  =  d>A;2wt2(LA  +  Aa)  -  d>A;2wt2LA  +  2^  ^^  -  Aa) 

+  2wt2(LC  "  Aa)  "  (20 c  iwtiLc  +  0C  2wt2Lc)  (11) 


AU  -  wjrf>A^t2  -  (0c?2t2  +  20C,ltl)JAa 
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or 


Aa-0  A^  '  HT  "  W[^)A,2t2  "  (^0,2*2  +  2d>C,ltl). 
Substitution  of  equations  (8)  and  (13)  into  equation  (1)  yields 

G  =  P(eA  "  ec)  "  w_^)A,2t2  "  (^0,2*2  +  2<^)C,ltl)_ 


This  equation  is  evaluated  in  the  appendix  in  terms  of  applied  stress,  temperature  change, 
material  parameters,  and  specimen  configuration  and  leads  to 


rc 

(2‘lEl  +  t2E2>r 

;  -  ATjfflj  -  Og) EgJ 

where  the  temperature  change  AT  is  the  difference  between  the  debond  test  tempera¬ 
ture  and  the  cure  temperature  of  the  adhesive. 

Correlating  Parameters 

In  the  previous  section  the  expression  for  the  strain -energy  release  rate  G  was 
derived  on  the  implicit  assumption  that  the  debond  advanced  while  the  specimen  was  sub¬ 
jected  to  a  constant  stress  condition  (constant  applied  stress  and  temperature).  However, 
for  cyclic  debonding  the  relation  between  stress  and  debond  extension  is  unknown.  In  fact, 
the  debond  may  or  may  not  extend  throughout  the  entire  stress  cycle.  Consequently,  the 
actual  value  of  G  for  cyclic  debonding  is  unknown.  However,  the  aforementioned  anal¬ 
ogy  between  cyclic  debonding  and  crack  propagation  in  metals  suggests  that  the  debonding 
process  may  be  correlated  with  the  maximum  cyclic  stress,  the  range  of  cyclic  stress, 
or  with  a  combination  of  the  two.  Consequently,  three  G  parameters,  each  incorporating 
one  of  these  three  stress  conditions,  are  developed  herein  and  evaluated  as  correlation 
parameters. 

To  develop  the  G  parameters  for  each  of  the  three  stress  conditions,  equation  (15) 
was  rewritten  in  the  form 

g  =  'mr.  *  He, 

E2(2tlEl  +  t2E2) 

where 


a  =  S  -  AT 


(“1  -  ^): 
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The  G  parameter  for  the  maximum  stress  condition  is  derived  by  substituting  amax 
for  a  in  equation  (16),  where 

amax  =  Smax  -  AT (a? j  -  Q'2)E2  (1^) 


Thus, 


rmax 


tit2ElW 


E2  2tiEi  +t2E2 V- 


Smax  -  AT  (o'!  -  a2)E2l 


(18) 


The  G  parameter  for  the  stress  range  condition  is  derived  by  substituting  Act  for  a 
in  equation  (16),  where 

Act  =  °rnax  ”  amin  “  Smax  -  ^min  (iO) 


and 


amin  ~  ^min  ~  AT(a?i  "  ^2)E2 


(20) 


Thus, 


tit2Eiw 


rano;e 


E2(2tlEl  +  t2E2) 


(Emax  -  in 


(21) 


The  G  parameter  for  the  combination  condition  is  derived  by  substituting  crmax  Act  for 
in  equation  (16),  where 


°max  Aa  ”  j^max  ~  ^(al  ~  Q?2)E2fSmax  _  ^min 


(22) 


Thus, 


t-itoEiw  r 

^max* range  jz — ~ - rf^max  -  AT(ai  -  Q?2)e2|  (®max  -  Smin)  (23) 

E2(2tlEl  +  t2E2)u 


These  three  parameters  are  evaluated  for  specific  materials  in  the  next  section. 
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RESULTS 


Cyclic  debonding  tests  were  conducted  for  each  of  three  material  systems  and  are 
discussed  in  the  following  sequence:  First,  typical  test  results  are  discussed  to  illustrate 
the  manner  in  which  cyclic  debonding  progressed.  Second,  the  topographies  of  the  failure 
surfaces  are  analyzed.  And  third,  the  degree  of  correlation  of  the  debonding  data  with 
each  of  the  three  strain -energy  release  parameters  developed  in  the  previous  section  is 
evaluated. 


Typical  Debonding  Behavior 

Figure  4  shows  a  typical  variation  of  debond  length  with  applied  load  cycles.  In  all 
cases  debonding  was  initially  nonlinear  but  became  linear  after  the  debond  progressed  a 
short  distance. 

For  all  three  material  systems,  nonlinear  debonding  was  confined  to  a  region  near 
the  change  of  cross  section  of  the  specimens.  This  nonlinearity  appears  to  be  related  to 
a  transition  of  failure  mode.  In  figure  4,  for  example,  the  light  portion  indicates  failure 
within  the  adhesive  and  the  darker  portion  indicates  failure  within  the  composite  matrix. 
(Herein,  matrix  refers  to  the  epoxy  of  the  composite  core.)  Early  transition  of  failure 
modes  occurred  in  all  three  material  systems. 

Because  nonlinear  debonding  occurred  only  during  a  small  portion  of  each  test,  the 
correlation  analysis  presented  herein  considers  only  the  significant  linear  portion. 

Examination  of  Failure  Surfaces 

The  possible  failure  modes  for  the  specimens  tested  are  numerous;  failure  can  occur 
within  the  adhesive,  at  the  interfaces  between  the  aluminum  and  the  composite,  or  in  the 
composite  itself  at  the  fiber-matrix  interfaces.  Because  these  modes  are  intimately  re¬ 
lated  to  the  debond-propagation  rates,  the  failed  surfaces  of  the  specimens  were  examined 
with  a  scanning  electron  microscope.  The  results  of  these  examinations  are  discussed  in 
the  following  paragraphs. 

Material  system  I  (fig.  5).-  When  this  system  failed,  adhesive  remained  on  both  the 
composite  and  the  aluminum  sheet.  Within  the  adhesive,  separation  occurred  at  the  inter¬ 
face  between  the  scrim  cloth  and  the  adhesive.  The  figure  shows  a  typical  failure  surface 
as  viewed  toward  the  debonded  aluminum  sheet.  At  the  lowest  magnification  (fig.  5(a)),  the 
failure  surface  appears  to  have  a  woven  pattern.  This  pattern  is  the  imprint  of  the  scrim 
cloth  on  the  adhesive.  (After  debonding  the  scrim  cloth  remained  with  the  composite 
core.)  Figure  5(b)  shows  these  imprints  at  a  higher  magnification.  At  the  highest  magni¬ 
fication  (fig.  5(c)),  the  individual  scrim  cloth  fiber  impressions  exhibit  striations.  In  this 
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system,  the  striations  were  difficult  to  find.  However,  striations  were  more  apparent 
in  the  other  material  systems  and  will  be  discussed  later. 

Material  system  II  (fig.  6).-  This  system  failed  primarily  at  the  fiber -matrix  inter¬ 
face,  with  some  failure  at  the  matrix-adhesive  interface.  After  failure,  the  adhesive 
remained  on  the  aluminum  sheets.  The  figure  shows  a  typical  failure  surface  as  viewed 
toward  the  debonded  aluminum  sheet.  At  the  lowest  magnification  (fig.  6(a)),  the  smooth 
areas  are  fiber -matrix  interface  failure  surfaces  and  the  coarse  areas  are  matrix- 
adhesive  failures.  On  closer  examination  (fig.  6(b)),  impressions  of  the  graphite  fibers 
can  be  seen  on  the  smooth  failure  surface.  At  a  higher  magnification  (fig.  6(c)),  a  graphite 
fiber  impression  exhibits  striations.  These  were  observed  for  the  entire  range  of  rates 
measured  with  this  material.  The  striation  spacings  were  in  general  agreement  with  the 
debonding  rates. 

Material  system  III  (fig.  7).-  In  contrast  to  material  system  II,  this  system  failed 
primarily  at  the  matrix -adhesive  interface,  with  some  failure  at  the  fiber -matrix  inter¬ 
face.  After  failure,  the  adhesive  remained  on  the  aluminum  sheet.  The  figure  shows  the 
failure  surface  as  viewed  toward  the  debonded  aluminum  sheet.  At  the  lowest  magnifica¬ 
tion  (fig.  7(a)),  the  coarse  area  is  the  matrix -adhesive  failure  surface  while  the  smooth 
area  is  the  fiber-matrix  failure  surface.  On  closer  examination  (fig.  7(b)),  S-glass  fiber 
impressions  can  be  seen  on  the  smooth  surface.  At  a  higher  magnification  (fig.  7(c)),  the 
S-glass  fiber  impressions  exhibit  striations.  As  in  material  system  II,  these  striations 
were  observed  in  all  tests  and  generally  correspond  to  the  measured  debonding  rates. 

Correlation  of  Debond-Propagation  Rates 

The  measured  debond-propagation  rates  da/dN  and  the  corresponding  Smax  val¬ 
ues,  thicknesses,  and  stress  ratios  are  shown  in  tables  II,  III,  and  IV.  For  each  material 
system,  at  least  two  specimen  thicknesses  were  tested  for  each  of  several  stress  ratios 
ranging  from  0.10  to  0.48.  In  the  following  discussion,  the  propagation  rates  for  each 
material  system  are  plotted  against  the  three  previously  developed  strain -energy  release 
rate  parameters  Gmax,  Grange,  and  Gmax.range  (eqs.  (18),  (21),  and  (23), 
respectively). 

Material  system  I.-  In  figure  8(a)  the  rate  data  are  plotted  against  Smax.  The  sym¬ 
bol  shapes  indicate  different  thicknesses  and  the  symbol  shadings  indicate  various  ranges 
of  stress  ratio.  The  layering  of  these  data  indicates  that  for  a  given  Smax,  the  thicker 
specimens  debonded  faster  than  the  thinner  ones,  and  the  specimens  tested  at  lower  stress 
ratios  debonded  faster  than  specimens  tested  at  higher  stress  ratios. 

Figures  8(b),  (c),  and  (d)  show  the  data  for  material  system  I  plotted  against  the 
three  G  parameters  developed  earlier.  Figure  8(b)  shows  the  data  plotted  against 
Gmax  (eQ-  (18));  figure  8(c)  shows  the  data  plotted  against  Grange  (eq.  (21));  and  fig- 
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ure  8(d)  shows  the  data  plotted  against  Gmax. range  (eq.  (23)).  In  all  these  figures,  an 
equation  of  the  form 

=  c(G)n  (24) 

dN 

was  fit  to  the  data  using  a  least -squares  fit.  (In  the  least -squares  method  the  arbitrary 
constants  c  and  n  are  chosen  to  give  the  smallest  error  between  values  from  eq.  (24) 
and  test  values.)  Mathematically,  the  goodness  of  fit  is  represented  as  the  sum  of  the 
squares  of  the  errors  2r2  between  the  calculated  values  and  the  test  values  -  the 
smaller  Sr2,  the  better  the  fit.  Visually,  the  goodness  of  fit  can  be  represented  by  con¬ 
fidence  limits  as  discussed  in  reference  5.  The  95-percent  confidence  limits  are  shown 
in  the  figures  by  solid  lines  -  the  closer  these  lines  lie  together,  the  better  the  correla¬ 
tion.  The  95-percent  confidence  limits  used  herein  indicate  a  95-percent  probability  that 
the  true  mean  of  the  test  data  lies  within  the  limits.  For  material  system  I,  Grange 
best  correlates  the  data. 

Material  system  II.-  Figure  9(a)  shows  the  debond -propagation  rates  for  two  speci¬ 
men  thicknesses  tested  at  several  stress  ratios  plotted  against  Smax.  Again,  for  a  given 
value  of  Smax?  lower  stress  ratios  and  thicker  specimens  produced  faster  debonding. 

Figures  9(b),  (c),  and  (d)  show  these  rates  plotted  against  Gmax,  Grange,  and 
^max-range*  As  f°r  material  system  I,  Grange  correlates  the  data  better  than  Gmax. 
But,  unlike  material  system  I,  Gmax.range  correlates  the  data  better  than  Grange 
alone.  In  this  case  da/dN  appears  to  be  a  function  of  both  maximum  load  and  load 
range. 

Material  system  HI.  -  Figure  10(a)  shows  the  debond-propagation  rates  for  material 
system  III  plotted  against  Smax.  As  with  material  systems  I  and  II,  lower  stress  ratios 
and  thicker  specimens  produced  faster  debonding. 

Figures  10(b),  (c),  and  (d)  show  the  rates  correlated  as  before.  The  data  are  corre¬ 
lated  best  by  Gmax.range,  and  da/dN  appears  to  be  a  function  of  both  maximum  stress 
and  stress  range.  But,  in  contrast  to  material  systems  I  and  II,  Gmax  correlated  the 
data  better  than  Grange. 

To  summarize,  for  all  three  material  systems  an  equation  of  the  general  form 
c  _  t1t2E1w  ^ 

^2(^1%  +  *2^2) 

(which  is  eq.  (16))  correlated  the  data  quite  well.  However,  the  specific  form  of  this  equa¬ 
tion  that  best  correlated  the  data  seemed  to  be  material  dependent.  For  material  system  I, 


11 


Grano-e  best  correlated  the  data;  but  for  material  systems  II  and  III,  Gmax-range  best 
correlated  the  data.  The  coefficients  and  exponents  determined  by  fitting  equation  (24)  to 
the  test  data,  and  the  sum  of  residuals  for  all  three  material  systems  are  shown  in  the 
following  table: 


Correlating  parameters 

Material  system 

i 

II 

III 

c 

n 

2r2 

c 

n 

Sr2 

c 

n 

Hr2 

Gmax 

5.94  ~14 

3.94 

16.7 

9.55"12 

1.81 

0.48 

2.61-15 

3.91 

1.62 

Grange 

5.88-13 

3.73 

5.59 

6.13 "10 

1.28 

.18 

6.31-12 

2.75 

2.18 

Gmax- range 

4.83"14 

4.22 

7.84 

4.56-H 

1.72 

.14 

1.17-13 

3.40 

1.09 

Because  the  number  of  data  points  varied  among  the  systems,  Sr2  for  different  systems 
should  not  be  compared.  The  finding  that  different  forms  of  G  correlated  the  different 
material  systems  may  be  attributed  to  a  number  of  factors:  the  assumed  forms  of  the 
correlating  parameters  and  of  equation  (24),  the  observed  differences  in  failure  mecha¬ 
nisms  among  the  material  systems,  the  assumed  value  of  AT  in  material  systems  II 
and  III,  and  so  forth.  Evidently  the  basic  key  parameters  and  functional  forms  are  yet  to 
be  defined.  Consequently,  variations  of  best  correlating  parameters  for  different  material 
systems  are  not  surprising,  and  empirical  correlating  expressions  applicable  for  one 
bonded  system  may  not  be  appropriate  for  another  system. 

CONCLUDING  REMARKS 

Cyclic  debonding  rates  were  measured  for  room -temperature  tests  of  three  mate¬ 
rial  systems:  graphite/epoxy  bonded  to  aluminum  at  room  temperature,  graphite/epoxy 
bonded  to  aluminum  at  an  elevated  temperature,  and  S-glass/epoxy  bonded  to  aluminum  at 
an  elevated  temperature.  Loading  in  all  tests  was  at  constant  amplitude.  During  each 
test,  debonding  was  monitored  with  a  photoelastic-coating  technique.  The  rates  of  debond¬ 
ing  were  compared  with  three  expressions  based  on  strain-energy  release  rate  in  terms 
of  the  maximum  stress,  stress  range,  and  a  combination  of  the  two. 

The  debond -propagation  rates  were  significantly  influenced  by  both  adherend  thick¬ 
ness  and  the  cyclic  stress  ratio.  For  a  given  maximum  stress,  the  thicker  specimens  de- 
bonded  faster  than  the  thinner  ones,  and  the  specimens  tested  at  lower  stress  ratios  de- 
bonded  faster  than  specimens  tested  at  higher  stress  ratios.  Microscopic  examinations 
of  the  debonded  surfaces  showed  significantly  different  failure  mechanisms  for  different 
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material  systems.  For  graphite  bonded  to  aluminum  with  a  room -temperature -curing 
adhesive,  failure  occurred  within  the  adhesive.  For  graphite  bonded  to  aluminum  with  an 
elevated -temperature -curing  adhesive,  failure  occurred  primarily  in  the  matrix  near  the 
adhesive  interface.  For  S -glass  bonded  to  aluminum  with  an  elevated-temperature -curing 
adhesive,  failure  occurred  primarily  at  the  composite -adhesive  interface. 

The  strain -energy  release  rate  correlated  the  data  quite  well  for  each  of  the  three 
material  systems,  but  the  form  of  the  best  correlating  expression  varied.  Consequently, 
empirical  correlating  expressions  applicable  to  one  bonded  system  may  not  be  appropri¬ 
ate  for  another  bonded  system,  and  caution  should  be  exercised  in  drawing  parallels  be¬ 
tween  different  material  systems  without  supporting  test  data. 

Langley  Research  Center 

National  Aeronautics  and  Space  Administration 
Hampton,  Va.  23665 
December  4,  1975 


13 


APPENDIX 


DEVELOPMENT  OF  STRAIN-ENERGY  RELEASE  RATE  EQUATION 

The  rate  at  which  strain  energy  is  released  when  a  metal  debonds  from  a  composite 
(fig.  3)  was  given  in  the  body  of  this  paper  (eq.  (14))  as 

G  =  P(eA  -  ec)  -  w  4>a, 2l2  ~  ($C,2i2  +  2d>C,ltl)  (A1) 

In  this  appendix,  the  right-hand  terms  are  evaluated  to  derive  an  explicit  expression 
for  G  (eq.  (15)).  The  strains  in  equation  (Al)  were  defined  by  requiring  equilibrium  of 
forces,  strain  compatibility  among  the  layers,  and  a  constitutive  relation  between  stress 
and  strain  in  specimen  regions  A  and  C.  Equilibrium  is  satisfied  by 

P  =  w  J  a  dy  (A2) 

For  the  one -dimensional  analysis  used  in  this  problem,  strain  compatibility  was  assured 
by  assuming  that  the  strain  through  the  specimen  thickness  was  constant.  The  constitutive 
relation  for  the  problem  is  given  as  (ref.  6,  p.  259) 

a  =  E(e  -  a  AT)  (A3) 

These  three  relationships  are  used  to  calculate  the  strains,  strain-energy  density,  and 
strain-energy  release  rate  in  regions  A  and  C,  as  shown  in  the  equations  that  follow. 

In  region  A,  equilibrium  is  satisfied  when 

P  =  st2w  =  o-A?2t2w  (A4) 

The  constitutive  equation  (A3)  for  region  A  is 

CTA,2  =  e2(6A  -  a2  AT)  (A5) 

Substituting  equation  (A5)  into  equation  (A4)  and  solving  for  eA  yield 

£A  =  A  +  Q!2at  (A6) 

e2 

In  region  C,  compatibility  is  satisfied  when 

eC,l  =  eC,2  =  eC  (A?) 
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Equilibrium  in  region  C  is  satisfied  when 

P  =  St2W  =  2t]Wcrc;;l  +  ^C,2^2W  (A8) 

The  constitutive  relationships  in  region  C  for  the  metal  and  composite,  respectively,  are 


^C,l  =  El(€C  -  al  AT) 

(A9) 

aC,2  =  E2(eC  “  a2  AT) 

(A10) 

Substituting  equations  (A9)  and  (A10)  into  equation  (A8)  and  solving  for 

lead  to 

St  2  +  ft  i^i  +  ^2^2q?2) 

C  2Eltl+E2t2 

(All) 

The  strain-energy  density  can  be  expressed  as 

,  _  a2  _  [E(e  -  a  AT)]2 

0  2E  2E 

(A12a) 

or 

<b  =  |[e2  -  2ea!  AT  +  a2(AT)2] 

(A12b) 

For  region  A,  substituting  equation  (A6)  into  equation  (A12b)  yields 


*A,2  “XL1 


6  A 


-  2e^O!2  AT  +  a 


:(AT)2] 


(A13) 


Similarly,  for  region  C,  substitution  of  equation  (All)  into  equation  (A12b)  yields  for  the 
metal 


6 


2 


2eCa  1 


AT  +  a 


(A14) 


and  for  the  composite 
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Substitution  of  equations  (A6),  (A8),  (All),  (A13),  (A14),  and  (A15)  into  equation  (Al) 

yields 


G  = 


tlt2Elw 


E2(2tlEl  +  t2E2, 


[S  -  AT(ai  -  a 2)^2] 


the  desired  expression  for  strain-energy  release  rate  (eq.  (15)). 


(A16) 
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TABLE  I.-  MATERIAL  PROPERTIES 


Material 

Modulus  of 
elasticity,  E, 
GPa 

Thermal  coefficient 
of  expansion,  a, 
jliK'1 

7075-T6  aluminum  alloy 

71 

22.50 

Graphite /epoxy 

131 

-.38 

S -glass/epoxy 

61 

_ 

3.60 

TABLE  II.-  RESULTS  OF  CYCLIC  DEBONDING  TESTS  ON  MATERIAL  SYSTEM  I 
(ALUMINUM  AND  GRAPHITE /EPOXY,  EA-927R  ADHESIVE  CURED  AT  RT) 


Thickness 

R 

Ciax* 

MPa 

(a) 

da/dN, 

(im/cycle 

(b) 

Cax’ 

j/m 

(c) 

Grange* 

J/m 

(d) 

Gmax- range* 
J/m 

(e) 

ll>  l2> 

mm  mm 

0.51  0.84 

0.10 

485 

0.000714 

15 

12 

14 

.79 

.10 

722 

.0920 

33 

27 

30 

.81 

.10 

804 

.146 

41 

33 

37 

1 

f  .81 

.10 

1206 

2.50 

93 

75 

83 

1.02  1.68 

.34 

431 

.00318 

24 

10 

16 

1.60 

.10 

438 

.0168 

24 

20 

22 

1.78 

.21 

468 

.0268 

29 

18 

23 

1.63 

.10 

479 

.0229 

29 

24 

26 

1.65 

.42 

493 

.00152 

31 

11 

18 

1.63 

.42 

495 

.00190 

31 

10 

18 

1.85 

.29 

518 

.127 

36 

18 

26 

1.85 

.48 

528 

.00927 

38 

10 

20 

1.63 

.10 

532 

.166 

36 

29 

32 

1.70 

.01 

572 

.335 

42 

42 

42 

1.73 

.19 

577 

.373 

44 

29 

35 

1.68 

.10 

578 

.373 

43 

35 

39  ' 

1.75 

.36 

586 

.0589 

45 

18 

29 

1.91 

.18 

656 

1.03 

58 

39 

48 

1.60 

.10 

6  57 

.498 

55 

44 

49 

1.68 

.10 

683 

1.32 

60 

49 

54 

l  1.60 

.10 

766 

1.78 

74 

60 

67 

1.60  2.41 

.10 

219 

.00689 

9 

8 

8 

2.51 

.10 

269 

.0315 

14 

12 

13 

2.51 

.10 

322 

.0569 

21 

17 

19 

2.34 

.42 

343 

.000678 

23 

8 

13 

2.49 

.10 

377 

.00785 

28 

23 

25 

2.52 

.24 

381 

.00424 

29 

17 

22 

2.49 

.10 

431 

.163 

37 

30 

33 

2.46 

.34 

447 

.0866 

39 

17 

26 

2.44 

.21 

495 

.254 

48 

30 

38 

2.46 

.42 

505 

.0290 

50 

8 

29 

2.54 

.10 

535 

1.32 

57 

46 

51 

2.39 

.29 

559 

.399 

61 

31 

43 

2.54 

.10 

589 

5.59 

69 

56 

62 

2.39 

.01 

596 

4.52 

69 

68 

69 

F  2.34 

.19 

614 

3.18 

73 

48 

59 

aStress  in  composite  in  region  A. 

Calculated  from  linear  portion  of  debond  length  vs  cycles  curve. 
Calculated  using  <7max. 

Calculated  using  amax  -  crmjn. 
e  Calculated  using  crmax  and  <7max  -  crm^n. 
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TABLE  IIL-  RESULTS  OF  CYCLIC  DEBONDING  TESTS  ON  MATERIAL  SYSTEM  II 
(ALUMINUM  AND  GRAPHITE /EPOXY,  AF-126  ADHESIVE  CURED  AT  394  K) 


Thickness 

da/dN 

jum/cycle 

h 

J 

t2’ 

R 

®max, 

MPa 

^max’ 

j/m 

Grange > 
J/m 

Gmax- range) 
J/m 

mm 

mm 

(a) 

■(b) 

(c) 

(d) 

(e) 

1.02 

1.65 

0.10 

371 

0.0150 

58 

14 

29 

1.65 

.10 

424 

.0193 

67 

19 

35 

1.57 

.10 

542 

.0302 

89 

30 

52 

1.63 

.30 

548 

.0226 

92 

19 

42 

1.57 

.30 

552 

.0294 

91 

19 

41 

1.57 

.30 

552 

.0306 

91 

19 

41 

1.63 

.10 

585 

.0635 

100 

35 

59 

i 

1.65 

.10 

638 

.0592 

113 

42 

69 

I 

1.63 

.10 

692 

.112 

126 

49 

79 

1 

1.57 

.10 

704 

.0899 

127 

51 

80 

1.60 

2.59 

.10 

318 

.0279 

77 

17 

36 

2.46 

.10 

434 

.0437 

106 

30 

57 

1 

1 

2.62 

.10 

475 

.0627 

122 

37 

67 

| 

2.51 

.10 

484 

.0650 

122 

38 

68 

2.54 

.10 

535 

.0734 

139 

46 

80 

2.39 

.42 

549 

.0254 

146 

21 

55 

2.62 

.42 

549 

.0508 

141 

20 

53 

2.54 

.10 

588 

.116 

158 

56 

94 

V 

2.46 

.10 

652 

.144 

178 

68 

110 

aStress  in  composite  in  region  A. 

^Calculated  from  linear  portion  of  debond  length  vs  cycles  curve. 
cCalculated  using  ~ 

^Calculated  using 
eCalculated  using 


°max- 
°max  "  °min- 


crn 


max 


and 


°max  ”  amin' 
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TABLE  IV.-  RESULTS  OF  CYCLIC  DEBONDING  TESTS  ON  MATERIAL  SYSTEM  III 
(ALUMINUM  AND  S- GLASS/E  POXY,  AF-126  ADHESIVE  CURED  AT  394  K) 


aStress  in  composite  in  region  A. 

^Calculated  from  linear  portion  of  debond  length  vs  cycles  curve. 
cCalculated  using  (7max. 

^Calculated  using  crmax  -  o-mjn. 

Calculated  using  amax  and  crmax  -  crmin. 


UNIDIRECTIONAL 
COMPOSITE  — 


S 


ADHESIVE 

ALUMINUM 


Figure  1.-  Specimen  configuration. 


'  ^NEW  DEBOND  FRONT 

(b)  After. 


Figure  3.-  Test  specimen  before  and  after  an  increment  of  debonding  Aa 
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26 


cured  at  RT.  Material  system  I 


Graphite  fiber  impression 
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cured  at  394  K.  Material  system  II. 


S-glass  fiber  impression 
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cured  at  394  K.  Material  system  III. 


m 


0.51 

r\ 

w 

1.02 

□ 

1.60 

A 

0 

<  R  <  0.15 

a 

0.15 

i  R  <  0.30 

■ 

▲ 

0.30 

<  R  <  0.50 

(a)  da/dN  against  S 


95% 

confidence 


(b)  da/dN  against  Gmax. 

Figure  8.-  Correlation  parameters  for  debond -propagation  rate  for  aluminum 
and  graphite/epoxy  specimen  with  EA-927R  adhesive  cured  at  RT.  Material 
system  I. 


rane:e  ’ 


(c)  da/ dN  against  ^x*ange  * 


(d)  da/dN  against  Gmax.range- 
Figure  8.-  Concluded. 


range  7  ' 

(c)  da/dN  against  Grange. 


■  ▲  0.30  2  R  <  0.50 

_ _ _ - _ - _ I _ - _ - _ 1 

50  150  250 

^max- range  ’ 

(d)  da/dN  against  Gmax.range. 

Figure  10.-  Concluded. 
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